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Mode coupling is studied in a parabolic index fiber with a lossy boundary 
and square cross section. Statistical deviations of the fiber axis from perfect 
straightness and random changes of its width are considered as causing 
mode coupling. The excess loss caused by these mode coupling mechanisms 
and the loss penalty incurred for a certain degree of narrowing of the 
impulse response are estimated. 

I. INTRODUCTION 

Multimode optical fibers whose cores have parabolic distributions 
of the refractive index, 1-3 



-«.(i-£a). 



(i) 



are of great practical interest for light transmission over long distances, 
since their delay distortion is much less serious than that of con- 
ventional clad fibers. 

Since no optical fiber can ever be produced free of random imper- 
fections, it is important to know how statistical irregularities of the 
fiber affect its performance. Random irregularities of the fiber axis 
and random changes of the effective width of the fiber cause coupling 
among its modes. The mode losses are functions of the mode number. 
Absorption losses tend to affect all modes in the same way. However, 
if we assume that the fiber boundary either consists of an absorptive 
material or is a rough surface that scatters light, we must expect that 
higher order modes, whose fields reach strongly into the neighborhood 
of the fiber boundary, suffer much higher losses than lower order 
modes that are confined to the vicinity of the fiber axis. Coupling of 
the low-order modes to the high-loss, high-order modes increases the 
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overall waveguide losses. One objective of our study of the effect of 
waveguide irregularities is thus the determination of the excess losses 
caused by mode coupling. 

The second objective of this study of waveguide irregularities con- 
sists in determining the impulse response of the fiber. In the absence 
of coupling, each mode transports a fraction of the total power at its 
characteristic group velocity. Since the group velocities of different 
mode groups are not identical, pulse distortion results. 2,3 Mode 
coupling has the beneficial effect of improving the impulse response of 
the fiber. It is thus of interest to determine how much reduction of 
multimode pulse distortion can be achieved by random bends and 
random width changes of the fiber. 

The effect of random bends on parabolic index fibers with circular 
cross section has been estimated in an earlier paper. 4 Pure diameter 
changes of a fiber with circular cross section leave modes with different 
circumferential symmetries uncoupled. Statistical irregularities are 
unlikely to result in pure diameter changes without distorting the 
circular fiber cross section. However, an analysis of more general 
distortions of a fiber with nominally circular cross section is difficult 
to perform. For this reason we discuss a fiber with parabolic index 
distribution (1) but with square cross section. It seems reasonable to 
expect that the performance of a fiber with square cross section is 
similar to that of a fiber with circular cross section. We expect to find 
the correct order of magnitude of the losses and impulse response of 
the round fiber by examining its close relative, the fiber with square 
cross section. In particular, it should be possible to assess the relative 
effect of random axis deformations as compared to random width 
changes. In a square fiber, changes of only one set of opposing walls 
leave groups of modes uncoupled from each other. This situation 
corresponds to the circular fiber with pure diameter changes that 
leave modes of different azimuthal symmetry uncoupled. By allowing 
both sets of opposing walls to change their separation randomly, we 
are sure that all modes are coupled to each other. This model corre- 
sponds to a nominally round fiber whose cross section is deformed in 
an arbitrary way that does not conserve the circular symmetry. 

II. THE MODES OF THE PERFECT FIBER WITH SQUARE CROSS SECTION 

The modes of an infinitely extended medium with the distribution 
n* = n%(l-2^£\ (2) 
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of the square of the refractive index have the form 8 

* M ~ (n ir2" + 'p!g!)*w; " w 

It is 

r 2 = x 2 + y 2 . (4) 

The parameter a is an arbitrary constant that, in conjunction with A, 
determines the transverse dependence of the refractive index distribu- 
tion. However, in the round fiber it is convenient to associate a with 
the radius of the fiber boundary so that A is the relative difference 
between the values of the refractive index on axis and at the fiber 
boundary. 

The square of the refractive index distribution (2) does not follow 
precisely from (1). However, if one equation is regarded as the precise 
distribution of the corresponding quantity, the other holds approxi- 
mately provided A is small and we limit r to the range r ^ a. H p and 
H, are Hermite polynomials of degree p and q, and P is the p owe r 
carried by the mode. The parameter w is defined as 5 (k = uVcomo) 

V2a V 

w 



"(wtTl) (5) 



and determines the radius of the field distribution with p = q = 0. At 
r = w the field has decayed to 1/e of its value on axis. E P9 represents 
the transverse component of the electric field vector. The longitudinal 
field components are relatively much weaker and are not being con- 
sidered. The field (3) is only an approximate solution of Maxwell's 
equations. The propagation constant of the mode is given as 6 

= ,*„ = no* [l-?g( P + «+!)]'. (6) 

The modes of the square-law medium are mutually orthogonal and 
satisfy the relation 



2k \mo J- 



E pgEfj/jdxdy = P8 pp -8 qQ -. (7) 



So far, the fiber boundary has been ignored. The mode field (3) is an 
(approximate) solution of the guided-wave problem if we assume that 
the distribution (2) extends to infinity. However, each mode decays 
very rapidly outside of a certain region. For a given value of p the 
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field oscillates as a function of x passing through p zero crossings. The 
shape of the function 



ff «(^s) 



g-.«/»» 



is shown in Fig. 1 . At the point (see appendix) 



= x = w\p-\- 



(8) 



(9) 



the oscillatory behavior of the function changes to a rapid decay. 
If x' < a the presence of the wall does not interfere appreciably with 
the field distribution. However, if x' > a the field distribution is 
severely altered by the presence of the wall. Since we are assuming 
that the interaction of the field with the wall causes power dissipation 
either by absorption or by radiation, we consider those modes whose 
fields reach the vicinity of the wall with high field intensity as being 
effectively cut off. By replacing x in (9) with a we obtain the condition 
for the maximum value of p that can be allowed for low-loss modes. 

*- = ($~\~ nM 4\-\' (10) 

Since we are assuming that the boundary of the fiber has a square cross 




Fig. 1 — Plot of the function given by eq. (8). 
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section, we also must impose the same "cutoff" condition on q, 

We use the modes (3) of the infinite square-law medium (2) to 
describe the modes of the fiber with square boundary if p < p e and 
q < q c . If either p > p e and/or q > q c we regard the modes as so lossy 
that they are effectively cut off. This procedure is an approximation, 
but it allows us to obtain estimates (whose errors are unknown) to a 
complicated problem. 

III. COUPLING COEFFICIENTS 

The coupling coefficients between two modes are defined by the 
general expression 6,7 

K m ,, t = sJj /" *z£ dy(n> - n«)2S,Xv. ( 12 ) 

A square-law fiber with random axis deformations can be described 
by the following distribution of the square of its refractive index. 



* 2 = ^l 1 " l [(a! " /)2 + {V " ff)2]S } 



(13) 



We assume that / and g are both random functions of z and that 
f/a and g/a are small quantities. The deflection of the optical axis of 
the square-law medium has a far more important effect on the modes 
than the corresponding deviation of the fiber boundary. The deflection 
of the fiber boundary that results from the random bends of its axis 
is neglected. 

Substitution of (2), (3), and (13) into (12) results in 

J^.-^VF+l/W (14) 

and 

k-.... ± . = =£* ViTT»M. (15) 



ia- 



All other coupling coefficients vanish. The two choices (1 or 0) that 
are indicated under the square-root sign in (14) and (15) belong to 
the corresponding upper or lower sign of the subscript on the left-hand 
side. Random deformations of the fiber axis couple only neighboring 
modes. 
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For random width changes of the fiber we use the distribution 

ft, " n3 ( 1 - 2 [ciT^ + (A?]4 (16) 



For small values of f/a and g/a we can write this expression approxi- 
mately as follows : 

n» = n§ (l-|[^(l-2/)+^(l-2f)]a|. (17) 

For random changes of the width of the guide we obtain from (2), 
(3), (12), and (17) 



and 



£».**.« - ^^ V(p±l)(p + »/(*) (18) 



n kw 2 A 



*«■*••*» = ^T ^(9±D(? + o 2 )^(2). (19) 



All other coupling coefficients vanish. There are nonvanishing diagonal 
elements in this case. However, diagonal elements couple each mode 
only to itself. This self-coupling is of no importance if j{z) and g(z) 
have Fourier spectra with no zero (spatial) frequency component. 

IV. COUPLED POWER THEORY 

Mode coupling in waveguides with random irregularities can be 
described by coupled power equations. 8 

HP 1 HP N 

-J-! + ,7 -^ = -aj>, + £ M*% - P.)- (20) 

oz v„ at M = i 

P, is the average power carried by the mode labeled v, v, is its group 
velocity, and a v its power loss coefficient. The mode label v is used as 
an abbreviation for the set of labels p, q. The power coupling coefficient 
h vil is defined as follows: 8 

h vii = \K VI1 \F(0 V - 0,). (21) 

The coefficient i£„ M is the factor of the function f(z) or g(z) appearing 
in eqs. (14), (15), (18), and (19). The spatial power spectrum of the 
function /(z) [or 0(2)] is defined as 



F(6) =z(\J Q L fW^tfy ( 22 ) 
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It is assumed that L — >e° in (22). The symbol ( ) indicates an en- 
semble average. 

Since the random processes considered here tend to couple each 
mode only to one of its neighbors on either side (in mode label space) 
the equation system (20) can be converted to a partial differential 
equation whose variables are not only the length coordinate z and 
time t but in addition the two mode labels p and g. 9il ° If the number of 
modes below the effective cutoff value is very large, the set of discrete 
modes can be regarded as a quasicontinuum. We write 

is h P q, p 'q'{r p 'g' "pq) 

p'.q' 

= hpq,p+bp,q{Pp+bp,q * pq) T ^P5. P— Ap,g \* p— Ap.O * Pit 

-\- hpq,p, q +&q{Pp,q+bq * pq) T "pa, p,q-bq\P p,q-&q "pq) 

*(^U h(p) %] + { * qy il hiq)a i} (23) 

The last step follows by considering the discrete mode labels as con- 
tinuous variables and replacing differences by differentials. The nota- 
tion h(p) and h(q) serves as a reminder that, according to (14) and 
(15), the coupling coefficients depend only on p if q is held fixed, and 
(18) and (19) show that they depend only on q if p is held fixed. We 
thus obtain the approximate partial differential equation 

The average mode power P is now regarded as a continuous function 
of z, t, p, and q. The group velocity v is a function of p and q. We have 
omitted the loss term. We consider the modes as lossless if the vari- 
ables p and q remain below the cutoff values (10) and (11) and as 
having infinitely high loss if cutoff is exceeded. This fact can be incor- 
porated into the theory as a boundary condition by requiring 

P(Pc q.) = 0. (25) 

It has been shown in Ref. 8 how the pulse propagation problem can be 
solved by means of a perturbation method if the solutions of (24) for 
the time-independent case are known. We thus consider the trial 
solution 

P(z, t, p, q) = U(p)V(q)e-°* (26) 
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and obtain by substitution into (24) 

<*Yh£[ m %] + (A »l&[™%] + '-°- (27) 

We separate this equation into two ordinary differential equations by 
introducing the separation constant k 2 : 

J[X»5] +7^17-0 (28) 



(Ap) 
and 



lh»S] 



+W-° <! *> 



V. CALCULATION OF THE STEADY-STATE POWER LOSS 

The equation system (28) and (29) together with the boundary 
condition (25) (and an additional one to be discussed later) defines 
an eigenvalue problem. The lowest order eigenvalue en has the 
physical meaning of the steady-state loss of the statistical power 
distribution. 8 This quantity is of interest since it determines the 
additional losses that are caused by the statistical irregularities of 
the fiber. 

For random deformations of the fiber axis we obtain the power 
coupling coefficient h(p) from (14) and (21). 

h(p) = K(a)p (30) 

with 

JC(0) = (25^2)V(0). (31) 

We assume that f(z) and g(z) have identical power spectra so that 
h(q) follows from h(p) by replacing p with q. According to (6) the 
difference of the propagation constants of adjacent modes can be 
approximated as 

-Jo a 

/W« ~ 0m = n «"V- (32) 

This approximation is independent of the mode numbers. This means 
that only one spatial frequency (or actually a very narrow range of 
spatial frequencies) of the power spectrum F (ft) is responsible for 
mode coupling. For random axis deformations we have 

Ap = Aq = 1 (33) 
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so that we must solve the differential equation 

dp [ P dp J + K(Q) 
Its solution is a Bessel function of zero order, 



[ P f] + ^P = 0. (34) 



uw = j °( 2 ww^y <35) 

The choice of the Bessel function instead of a Neumann function, 
that would also solve (34), is dictated by an additional boundary 
condition. Since the partial differential equation (24) can be regarded 
as a diffusion process, we must require that no power diffuses into the 
lowest order mode p = from negative values of p. This requirement 
means that dP/dp = at p = 0. The solution (35) satisfies this 
condition. The solution of (29) is similarly 



v ^' J '( 2 tww^)- (36) 

The boundary condition (25) leads to 

VjT e = u, (37) 



Vtftf2) 
and 



sS^ =u - (38) 

The roots w, and u» are defined as solutions of the equation 

Jo(ttr) = 0. (39) 

Since the eigenvalues depend on the labels v and m, we attach these 
labels to a and obtain from (37) and (38) (note, p c = q c ) 

The steady-state power loss, the lowest order eigenvalue an, follows 
from (5), (10) (neglecting the term |), (31), and Ui = 2.405 

axis deformation: «r u = 5.78 -. F(O). (41) 

a* 

We have thus rederived the loss formula (42) of Ref . 4. 
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For random diameter changes we obtain from (18), (21), and (31), 
considering that the spacing (in /9-space) between adjacent coupled 
modes is now twice as large, 

h(p)=l^ 2 K(2Q)p\ (42) 

Since the number of modes is assumed to be large, we have used the 
approximation p(p — 1) tt p 2 . With Ap = Aq = 2 we obtain from 
(28) and (42) 

The solution of this differential equation is 

U = -j= cos (p. In p + <t>v) (44) 

Vp 



with 



- /*_*__ 1. (45) 



The solution of (29) is correspondingly 



1 

7 = ^= cos (p. In q + <P M ) (46) 



with 

/a 2 g - k 2 1 
p " ylw*K(2tt) 4 



(47) 



These solutions have a singular behavior at p = or q = 0. However, 
we must keep in mind that p and q are really discrete quantities. 
Considering them as continuous variables is an approximate procedure 
that can work only for very large values of p or q where the relative 
difference between adjacent discrete values becomes small. Since 
In p = for p = 1, we allow p and q to vary only between 1 and p c . 
The requirement that no power diffuses across the lower limit of the 
range of the variables imposes the conditions 

^ dp y,_i " 

and (48) 

\ ~dq A=i ~ 

These conditions lead to the determination of the phase terms via 
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the equations 

tan <p, = - h- ( 49 ) 

and 

tan</> M =- t^-- (50) 

2p„ 

The boundary condition (25) leads to 

p, = ,— [ (2" - 1) £ + arctan i- 1 • (51) 

In p c L " ^P" J 

p„ as well as p M are solutions of this equation with integer values of v. 
Since the values of p, are now known, we obtain the eigenvalue o th 
from (45) and (47) 

^ = ^K(2n)(p 2 , + p* + ^). (52) 

For the lowest order eigenvalue, that is, for the steady-state power 
loss coefficient, we obtain with the help of (5) and (31) 

width changes: cm = (4p? + 1) ^ f (20). (53) 

The solution of (51) is not a constant. It depends on the waveguide 
parameters through its dependence on p c = (a/w) 2 . A plot of pi as a 
function of p e is shown in Fig. 2. 

The forms of the steady-state power loss coefficients (41) and (53) 
are very similar. The power spectra describe the deflection of the fiber 
axis from its nominally straight position or the changes of the width 
of one of the transverse fiber dimensions. However, the excess loss 
caused by random changes of the width of the fiber depends on a 
component of the power spectrum at twice the spatial frequency com- 
pared to the excess loss for random bends of the fiber axis. 

VI. CALCULATION OF THE PULSE WIDTH 

The width of the impulse response of a multimode fiber that is 
long enough for the steady-state distribution to establish itself is 
given by the formula: 8 



» (Gu, VG V)t y 



M = 4 \ L E 



(54) 



The term with «r,„ = on is excluded from the sum. The functions 
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Fig. 2 — Plot of the parameter pi as a function of p e . 



G,. are defined as 



G„ = A,JJMVM 



(55) 



and V is by definition the difference between the inverse group velocity 
of the modes minus the inverse of the maximum group velocity. 4 



V = 



1 



c cjiok 2 a? 



(p + g) 2 - 



»(p, q) 

The expression in parenthesis is an abbreviated way of writing 



(0U, V0 W „) = fdpjdqGuVG,,. 



(56) 



(57) 



The integrals extend over the entire range of p and g variables from 
either or 1 to the cutoff value p c = q c . Requiring the normalization 



JdpjdqGl= 1, 
we have for random axis deformations : 

. Jo [ Vn>\\— ) «/o ( U?\ — ) 

G = w \ V pe ' V vQ</ 
'* a 2 JiMJiM 
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(58) 



(59) 



and for random width changes : 
A 



G„ u = ^ cos [p„ In p + </>,] cos [p„ In g + 0„] 



with 



A.. = 2 | [in p. + r^]['n«« + pjng] 



H 



(60) 



(61) 



The integrals (57) have the following solutions. For axis deforma- 
tions (p. 5^ 1) : 

(flu, VGi,) = (On, VG^i) 



32Ap 2 c uiw, 



c/ioWfr? - u2) 



2\2 



2u\ - 1 6(ti? + ul) 1 



3u? 



(i*i - ujj) 



and for v, p. ?± I 

{Gn, VQ,J - ^^- 2 (u? _ W 2 )2 (w? _ M 2 )2 Vc 

For random width changes (p ^ 1) : 



(62) 
(63) 



(Gu, VGi„) = (Gn, VG M i) =- 4 



cnokras 



PiPm 



pMnA 



u 



K-i) 



2p 

X ' (1 + 4p?)[l 4- (pi + Pm) 2 X1 + (oi - p,y] 
1 + 2(p? + pD 1 



<- 



1)" + 



pW(l+4p?)(l+4p 2 ) 

1 



A„[4 + (pi + p M ) 2 ][4 + (pi - p„) 2 ] 

5 + 2(p? + pi) 



X 



<- 



1)" + 



PW(1 +4p?)(l + 4p 2 ) 



1! 



(64) 



and for i>, p ?* 1 
(G 1X , VQ,J = 8 



cnok 2 a 2 



p'iAuA, 



PlPxPji 



X [1 + ( Pl + p,) 2 ][l + (pi " P,) 2 ][l + (pi + P„) 2 ][l + (Pi - Pm) 2 J 
1 + 2(p? + pl) 



X 



(-1)' + 



PeV(l +4p?)(l +4 P 2 ) 



X 



(-1)" + 



(65) 



1 + 2(p? + p 2 ) 
PcV(l+4p?)(l+4p 2 ) 
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Evaluation of (54) with the help of (62) and (63) yields, for random 
deformations of the fiber axis, 



At = 



0.42Ap 2 e VL a 



(66) 



Equation (66) determines the pulse width of an impulse after it has 
traveled a distance L (L must be large enough so that the pulse has 
settled down to steady state) in the presence of random deformations 
of the fiber axis. The pulse width for uncoupled modes is obtained 
from (56) 

L L 4AL , 



AT = 



v (Pc, Qe) v(0, 0) cnok?d- 



Pc 



(67) 



The relative improvement of the width of the impulse response caused 
by mode coupling is characterized by the ratio 8 



R = 



At 



0.105 a 



&T VLtf(a) w 



(68) 



Mode coupling not only shortens the width of the impulse response, 
but it also leads to excess loss. In order to find out how much excess 
loss is associated with a given improvement of the width of the impulse 
response, we form the product of (41) with the square of (68) 



RtffnL = 0.032. 



(69) 
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Fig. 3— Plot of the function H(p e ). 
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Fig. 4 — Plot of the loss penalty as a function of p e . 

For random axis deformations, the product of the square of the im- 
provement factor with the excess loss is independent of the waveguide 
parameters and the statistics of the random axis deformations. 
For random width changes we obtain similarly 

The function H(p c ) is plotted in Fig. 3. The improvement factor is 

R= H^_a_ (71) 

VLlf(2a) w 

Finally, we obtain the loss penalty from 

R*a n L = (2 P 5 + 0.5)# 2 (P«)- (72) 

This function is shown graphically in Fig. 4. 

VII. DISCUSSION AND NUMERICAL RESULTS 

We have derived expressions for the steady-state loss and the loss 
penalty of graded index fibers with square cross section for the case of 
random axis deformation and random width changes. The most 
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conspicuous difference between these two types of fiber imperfections 
is the fact that, whereas the Fourier components of the function f(z) 
(describing the fiber axis) at the spatial frequency 12 are instrumental 
in the mode mixing process, the Fourier components at twice the 
spatial frequency, 212, determine the mode mixing process in case of 
random changes of the width of the fiber. This behavior can easily be 
understood. Consider a Gaussian beam of arbitrary width that is 
injected into the fiber off axis. 11 The beam undulates periodically 
around the fiber axis and also changes its width periodically. The 
undulations around the fiber axis have a period 11 

A = ^S? (73) 

while the width changes repeat themselves with half that period or at 
twice the spatial frequency. 11 Random displacements of the fiber axis 
couple to the deflections of the beam from its on-axis position. This 
deflection is driven by a Fourier component at the spatial frequency 

= £ = ^. (74) 

A a 

The beam width changes are correspondingly driven by changes in 
the gradient (width) of the parabolic index medium. It is thus clear 
that they respond to twice the spatial frequency. 

In order to be able to associate specific rms deviations of the fiber 
axis or rms width changes with fiber loss we have to consider a par- 
ticular statistical model. We choose (arbitrarily) a Gaussian correlation 
function 

</(«)/(* + u)> = oH-^\ (75) 

a is the rms deviation of the function f(z) and D its correlation length. 
The power spectrum of f(z) is known to be 12 

F(8) = Vtfif ^ I>e-<« , / ^ >^ (76) 

For a given value of 6 this function assumes its maximum value 

rcvM-i V2^e-°- 6 1.52 . . 

[_F(8)] D=Dm = ^ a 2 = -^- <r 2 (77) 

at 

An = V2/0. (78) 

Let us consider a numerical example. We use the following fiber 
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parameters 

a = 4.85 X 10- 3 crrn 

no =1.56 \- (79) 

A = 0.014 J 

At X = 1 Mm wavelength we have a/w = 6.3 or p e = 39.7. The 
difference between the propagation constant of adjacent modes is, 
according to (32), ft = 34.5 cm -1 . With = ft we calculate the excess 
loss values at the peak of the power spectrum at the value of the 
correlation length given by (78), D m = 0.041 cm. For random devia- 
tions of the waveguide axis we obtain from (41), (77), and (79) 
(a in cm) 

an = 6.44 X 10 6 ^ 2 (cm" 1 ). (80) 

In order to keep the excess loss below 10 dB/km = 2.3 10" B cm" 1 we 
must keep the rms deviation of the fiber axis below j = 2X 10~ 6 cm. 
However, this stringent tolerance requirement results from our as- 
sumption that the correlation length of the random irregularities of 
the fiber axis assumes its worst possible value (78). If, for example, 
the correlation length happens to be D = 0.5 cm we obtain instead 
of (80) 

an = 6.4 X 10" 2B a 2 (cm- 1 ) (81) 

SO that we can now tolerate 9 = 6 X 10 9 cm in order to keep the excess 
loss below 10 dB/km. This example shows that it is impossible to 
predict the excess loss to be expected from a practical square-law fiber 
unless the statistics of its irregularities are known precisely. 

For reasons of comparison we state the corresponding value for 
random width changes. In this case the spatial frequency that is 
instrumental in providing mode coupling is 2ft = 69 cm" 1 . The worst 
possible correlation length is now D m = 0.02 cm. From (53) and Fig. 2 
with p c = 40 we obtain (a in cm) 

an = 1.39 X 10 6 <f 2 (cm- 1 ). (82) 

The tolerance requirements of random width changes appear a little 
less stringent than those of random axis deformations. However, we 
have already pointed out that the excess loss value depends critically 
on the actual statistics of the fiber. Since the excess loss caused by 
random axis deviations depends on a different spatial frequency than 
the excess loss caused by random width changes, a loss comparison 
of the two effects is not possible. 
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Next we discuss the loss penalty that is incurred for a given im- 
provement of the width of the impulse response of coupled mode 
operation compared to uncoupled mode operation. The equations (69) 
and (72) show that the loss penalty is independent of the statistics of 
the fiber irregularities. This feature makes the loss penalty a useful 
quantity. In case of random variations of the fiber axis, the loss penalty 
is even independent of the fiber parameters and is simply a dimension- 
less number. Let us assume that we want to achieve a ten-fold improve- 
ment of the width of the impulse response compared to the impulse 
response of uncoupled multimode operation. In this case we have 
R = 0.1 and obtain from (69) for random deviations of the fiber axis 

eruL = 3.2 = 14 dB. (83) 

The length L needed to incur this loss and at the same time to achieve 
R = 0.1 depends on the statistics of the irregularities. However, eq. 
(83) tells us that it costs 14 dB in excess loss to achieve a ten-fold 
relative pulse width improvement. 

For random width changes, the situation is slightly different. Here 
the loss penalty depends somewhat on the fiber parameters. For the 
values used earlier we find from Fig. 4 with p c = 40, 

BVuL = 10" 2 . (84) 

The loss penalty for R = 0.1 is more favorable in this case, 
ail L = 4.34 dB. 

Fiber irregularities can be introduced intentionally in order to 
improve the impulse response. In the conventional fiber with a round 
core of constant refractive index that is surrounded by a cladding 
with constant index, the loss penalty for pulse distortion improvement 
can be reduced (in principle avoided) by tailoring the shape of the 
power spectrum carefully. 8 The reason that the shape of the power 
spectrum has an influence on the loss penalty is explained by the 
observation that the spacing (in /3-space) between adjacent modes of 
the conventional fiber is dependent on the mode number, so that a 
band of spatial frequencies of the power spectrum takes part in the 
mode coupling process. 

In case of the parabolic index fiber, only one spatial frequency 
(or actually a narrow range of spatial frequencies) is responsible for 
mode coupling. The shape of the power spectrum is thus immaterial, 
only its value at the spatial frequency 12 enters into the picture. The 
expressions (69) and (72) show that the loss penalty of the parabolic 
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index fiber is independent of the power spectrum. No loss advantage 
is to be gained by using especially shaped power spectra. One might 
think that an advantage could be gained by departing from the 
square-law index profile in order to change the mode spacing and 
sample more of the power spectrum. But as soon as the index distribu- 
tion deviates slightly from the parabolic profile the uncoupled impulse 
response becomes much broader. The mode coupling mechanism would 
now have to work against a far less favorable (uncoupled) impulse re- 
sponse so that it seems unlikely that an advantage can be gained 

in this way. 

Finally, we consider an example of pulse width reduction by random 
irregularities. We can introduce intentional deviations of the fiber 
axis from perfect straightness in order to cause mode coupling. Since 
the coupling process must be random, we could use deformation func- 
tions f(z) and g(z) that are sinusoidal in shape but have a random 
phase. This introduces a power spectrum centered around the spatial 
frequency of the sinusoidal process having a finite width. Instead of 
pursuing this idea further, we assume that we have somehow created 
an axis deformation whose power spectrum reaches beyond the fre- 
quency J2 of (32). For simplicity, and to have a definite case in mind, 
we choose 



F(fi) = - 



*jL \e\ <2n 

2fi ' ' 



o |0| > 2n. 



(85) 



This power spectrum is flat from zero spatial frequencies to a cutoff 
value of = 2fi and zero for > 20. The rms deviation a of the fiber 
axis from a straight line appears in (85). Using the fiber parameters 
(79) we obtain from (68) (ff, L in cm) 

_ 6.9 X 10- . m 

A ten-fold improvement of the width of the impulse response 
(compared to the uncoupled case), R = 0.1, over a length of 
L = 1 km = 10 5 cm requires an rms deviation of the fiber axis of 
a = 2.2 X 10 -6 cm. We already know that we pay for this improve- 
ment of the impulse response with a loss penalty of 14 dB. Very slight 
random deviations from perfect straightness are already very effective 
in providing mode coupling and improving the width of the impulse 
response. 
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For random width changes we have to allow for a wider power 
spectrum. Letting the power spectrum again extend twice as far as 
the effective spatial frequency, 2S2 in this case, forces us to divide 
(85) by 2. We thus find from Fig. 3 and (71) 

R = 837 >< 10 " - (87) 

R = 0.1 and L = 10 5 cm requires d = 2.6 X 10~ 6 cm. 

APPENDIX 

The function 

H p (j2?\er**i»'eru** (88) 

describes the modes of a square-law medium defined by 

n(x) = n (l - ^ a) . (89) 

The associated ray problem can be described by a paraxial Hamiltonian 
of the form 13 

"-■&- •«■ (90) 

The quantum mechanical treatment of this problem leads to an ex- 
pression for the "energy" E of the ray that has the form 14 

E=-^- (91) 

We define the "turning point" of the light rays associated with the 
wave field (88) by the condition that p x , which is proportional to the 
slope of the light ray, must vanish. That means that the ray trajectory 
is tangential to the optical axis as the rays turn back in their path 
leading them away from the axis. Using p x = and equating (90) 
and (91) we find the following condition for the turning point: 

n{x) = |*. (92) 

The propagation constant of this two-dimensional mode field is 6 

'-—* !>-"§(> +i)]' (93) 
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Substitution of (89) and (93) into (92) leads with the help of (5) to 
the formula (9) for the turning point. 

The physical argument advanced here serves the purpose of denning 
the range in which the Hermite polynomial has an oscillatory behavior. 
This range is given by 

-x ^ x <, x. (94) 

Outside of this range the Hermite polynomial grows monotonically to 
infinite values. However, since the Hermite polynomial enters the 
mode field (88) only as a product with a Gaussian function, the mode 
field decays rapidly without oscillation outside of the range given 
by (94). 
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